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Abstract 

We discuss expectations for the total and inelastic cross-sections at LHC CM energies ^/s = 
7 TeV and 14 TeV obtained in an eikonal minijet model augmented by soft gluon /c^-resummation, 
which we describe in some detail. We present a band of predictions which encompass recent LHC 
data and suggest that the inelastic cross-section described by two channel eikonal models include 
only uncorrelated processes. We show that this interpretation of the model is supported by the 
LHC data . 
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I. INTRODUCTION 



Models for the high energy behaviour of the total cross-section hold a large interest since 
Heisenberg's early attempts to describe the production of mesons through a shock wave 
mechanism [1]. The recent measurement of the inelastic cross-section at cm. energy of 7 
TeV by the ATLAS p] collaboration, together with preliminary CMS results ^ and the 
expected results from other LHC groups, have renewed the interest in the subject [1H7]. 
In addition, very recently, the TOTEM collaboration has released values for the total, the 
elastic and the inelastic cross-sections ^ and the Auger experiment has given a measurement 
of the p — air cross-section at y/s = 57 TeV 

The interest arises because total cross-sections access the large distance regime of particle 
interactions, i.e. a region where a QCD description is still lacking. In this paper, we present 
estimates for the total and the inelastic cross-sections at the LHC energies, ^/s = 7 and 
14 TeV, based on an eikonal minijet model augmented by soft-gluon resummation which we 
have developed over the years [TUHT2]. Our present aim is not to give newer fits to the data, 
but to use the model to further develop an understanding of the underlying large distance 
dynamics. This has led us to an interesting result, discussed in this paper, namely that the 
expression for the inelastic cross-section obtained from simple two channel eikonal models 
only describes the contribution of uncorrelated scattering processes. 

The eikonal model for high energy scattering [121 El provides a framework which in- 
corporates unitarity in the description of scattering processes and can be used to describe 
both the total and the inelastic cross-section. Following up on the earliest discussions of the 
role of low pt jets (aka mini-jets) in affecting the energy dependence of total cross-sections 
[T5HT7] . we have used the eikonal model to study the role of mini-jets in a QCD based frame- 
work, using as much as possible standard QCD phenomenological inputs, and techniques. 
Our model gives a description of the data on total cross-sections in terms of various QCD 
inputs at Leading Order (LO). It is able to address two major characteristics of the total 
cross-section, namely the sudden, power-like rise with energy at the ISR and the subsequent 
levelling off as needed by the Froissart bound. 

Our LO description of the rise is based on LO parton density functions and LO parton- 
parton cross-sections, implemented by soft gluon resummation of the initial state radiation, 
which in turn is parametrized through the same LO PDF used for the mini-jet calculation. 
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Non-leading effects in the resummation are incorporated in a non-perturbative ansatz for the 
singular, but integrable, coupling of soft gluons in the infrared region. Our model uses these 
inputs to obtain results for the total and the inelastic cross-section which are compared to 
the Tevatron and LHC energies. We notice that the inelastic cross-section poses yet another 
challenge, as the theoretical framework in which to model this quantity is not as easily 
defined as is the case for the total cross-section. We discuss this point in connection with 
eikonal models. 

The plan of this paper is as follows. In Sec. [IT] we discuss some features of the basic 
eikonal mini-jet model, and discuss the limitations of the two channel eikonal model. In 



Sec. |III| we briefly summarise the relevant details of our model, which includes soft gluon 
/ct-resummation, and in Sect. |I V| we present numerical results on the total cross-section and 
some of the quantities which crucially control its energy dependence. In Sect. |V] we apply 
our model to calculate the inelastic cross-section and discuss the range of model parameters 
required to obtain an adequate description of currently available data for both the total and 
inelastic cross-sections at 7 TeV and beyond. Our results are then compared with existing 



data and recently measured total and inelastic cross-sections at LHC. We end in Sect. VI 



with a comment on the limitations of eikonal models in describing total and inelastic cross- 
sections. 



II. TOTAL AND INELASTIC CROSS-SECTIONS IN EIKONAL MINI-JET MOD- 
ELS 

Eikonal mini-jet models are an extension of Glauber's theory of diffractive scattering 
[H], with the high energy dependence driven by QCD processess, such as production of 
low-pt(niini) jets. In addition to proton-proton scattering [181 IlH] and p — air processes 
[201 |2T] . the eikonal minijet models have been applied to photon processes [22l 123] . QCD 
inspired versions [211 |25] have been used to provide a unified description of proton and 
photon processes. Our version of the eikonal mini-jet model uses LO parton densities, 
determined from Deep Inelastic Scattering (DIS) and QCD evolution, and includes soft 
gluon fc^-resummation as a crucial component to access the very large distances probed by 
total cross-sections. Our model has been used to analyse total cross-sections for pp and 
pp [261 EZj, TriV and vrvr p8], [29] and 77 [301 EI], and is shown to give a satisfactory 
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description of the high energy behaviour of available data. 

Eikonal models are formulated in the b-space representation of the scattering amplitude 

F{s, t) = J d%f{b, s)=i j (fhe'^'-^l - e'^('''^)] (1) 

One then gets in impact parameter space 

d <7 elastic , i^((h,s)\2 /0\ 

Use of the optical theorem gives 

^ 2 / <?b*e[l - e'x.-.] . 2 / .=b|l - cos .)e--<M.l (3) 

With the above equations, and defining aimi = (^totai — (^elastic one obtains the total inelastic 
cross-section as 

cr,„,, = I d^l - e-2^™>^(M)] (4) 

The above equation presents the possibility of calculating the imaginary part of the eikonal 
function s), by relating it to the average number of independent inelastic collisions at 
given values of impact parameter b and cm. energy ^/s. It is in fact possible to obtain 
the above expression for the inelastic total cross-section through a semi-classical argument, 
based on the hypothesis that the scattering between hadrons takes place through multiple 
parton-parton collisions which are independently distributed. This corresponds to assuming 
a Poisson distribution around an average number of collisions n, namely 

P({n,n}) = i^ (5) 

One then calculates the scattering at each impact parameter value b between the scat- 
tering hadrons and integrates over the impact parameter space. Explicitly, 

a^nel{s) = [ ^(^^' = / ^'^[1 - e'^^^'^)] (6) 

n=l ^ 



Comparing Eq.(|4]) with Eq.(|6]) allows the identification n{b,s) = 2'^mx{b, s). Building a 
model for the average number of inelastic collisions then leads to Qmx and to making a 
model for the total cross-section as well. In particular, by putting 3?ex(&, s) ~ in Eq. (|3]), 
one gets 



and 

Such simple mini-jet models, as described above, can be used successfully to describe 
the total cross-section but the parameters, which allow to reproduce correctly the energy 
dependence of the total cross-section, predict inelastic cross-sections which are generally too 
small. Correspondingly, they predict elastic cross-sections too large and one needs further 
input to reproduce the differential elastic cross-section. An explanation, put forward in 
|32j . is that a simple description of the scattering in terms of just two channels, elastic 
and inelastic, or elastic and total, is not adequate and that (single and double) diffractive 
channels need to be treated separately. To avoid introducing more parameters, one can insist 
on a two-channel model, as in Eqs. ([6]), ^ and ([s]), but in such case the inelastic cross- 
section, calculated through these equations in the eikonal mini-jet model, would not include 
correlated processes. The reason for this comes from comparing Eqs. Q and (|6|, which 
show that the inelastic formula, Eq. (|4]), includes only those processes where the underlying 
dynamics contains no correlations between collisions. 

The above observation deserves a brief consideration. Many models, which use the eikonal 
framework, apply Eqs. Q and ^ to calculate the inelastic cross-section. The above dis- 
cussion shows that, in so doing, the result will miss the contribution of correlated processes, 
such as diffraction. This is particularly important when extracting total cross-section data 
from cosmic ray experiments, which measure the inelastic pp cross-section [33]. Another 
case, where the inelastic cross-section matters, is the calculation of the Survival Probability 
for Large Rapidity Gaps [3l]. Such models, [5l [35j or our own [26], use the probability for 
occurence of no inelastic collisions in 6-space, and caution must be exerted in interpreting 
what type of processes are excluded. 

We shall return to this point at the end of the paper and investigate it in detail in 
future work. In the meanwhile, we present the results we obtain with the model described 
by Eqs. ([6]), ([T]) and ^ and study the role played by the parameters in describing total 
and inelastic cross-sections. We shall then compare our results with the recent ATLAS 
Collaboration [2] analysis of inelastic collisions at LHC cm. energy ^/s = 7 TeV, as well 
as with the preliminary CMS results [3] and with the TOTEM value for the total cross- 
section [S]. 
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III. THE EIKONAL MINUET MODEL WITH SOFT /fef-RESUMMATION 



The average number n{b, s) of collisions can be written as the number of scattering 
centers per unit area A{b,s) multiplied by the scattering cross-section cr(s), schematically 
A{b, s)cr(s). In the parton model, we need both (i) the various parton-parton cross-sections 
as well as (ii) the parton densities for an evaluation of A{b, s), a{s) and hence n{b, s). 

In minijet models, the rise of the total cross section with energy is related to the increasing 
probability of perturbative low-x parton-parton collisions. Defining as Ptmin the threshold 
for perturbative treatment of parton-parton scattering, this partonic contribution to cr(s) 
can be calculated as 

dpt 



PtTl 



^i^ti^'^ Ptmin) = I dpt dXi dX2 22fi\AiXl,Pt)fj\B{x2,Pt) • (9) 

J4p2/s J4p2/(a;is) ■ ■ ^ 



i,j,k,l 

Here, i, j, k, I denote the partons and Xi,X2 the fractions of the parent particle momentum 
carried by the parton. = y/xiX2S and a are the center of mass energy of the two parton 
system and the hard parton scattering cross-section respectively. The numerical evaluation 
of this quantity strongly depends upon ptmin and the chosen set of Parton Density Functions 
(PDFs). 

The minijet model is moot regarding the region pt < Ptmin) and the underlying formalism 
is well defined only for perturbative values of momenta of the scattered partons for pt > Ptmin- 

rp J 1 . J J r 1 • 1 outqoinq parton , Tj. — • x 

io take mto account processes tor which Pt < Ptmin-, one can split n into a soft 

and a hard part. i.e. 

n(6, S) = Usoftib, S) + Uhardib, s) (10) 

We then factorize each term into an overlap function A{b, s) and a cross-section cr(s), as 

nsoft/hardib, s) = Asoft/hardib, s)asoft/hardis) (11) 

The subscript hard is meant to indicate the perturbative origin of the mini-jet contribution 
to <y totalis) in calculating n(6, s), as given by Eq. (|9]). To proceed further in the calculation of 
n(6, s), one needs to supplement this simple model with the knowledge of density of partons 
in impact parameter space as well as a model for figoftip, s). These issues will be discussed 
below and in the next section. 
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A. Overlap functions in the transverse space 



In the earliest discussions of eikonal models (TH], A{b, s) was taken to be independent of 
s and was obtained by taking the convolution of the Fourier transform of the form factors 
(FF) of the colliding particles A, B, given by 

A^Fib) = I (12) 

However, these simple parametrisations of the overlap function A{b,s), augmented with 
a model for nsoft{b,s) failed to reproduce correctly the observed energy dependence [TT] . 
When using in ajet standard parton densities determined from Deep Inelastic Scattering 
(DIS), one found an energy rise that was either too fast or started too early. An acceptable 
description of the data, from the beginning of the rise to the TeVatron energies, was not 
possible without modifying the parton densities in an ad hoc manner. 

We were able to obtain a good description of the energy behaviour using soft-gluon kt- 
resummation, in a model (BN) labelled after Bloch and Nordsieck [36]. It reflects the idea 
that any description of the scattering process between charged particles needs to include 
an infinite number of soft quanta in order to obtain a finite result. This model provides 
an energy dependent form factor, with the energy dependence introduced through the kine- 
matics of single soft gluon emission in the 2 — )■ 2 parton-parton collision. In the BN model, 
the parton distribution in 6-space, is determined as the Fourier transform of the probability, 
c/^P(Kt), that the initially coUinear parton-parton pair acquire, in the collision, an overall 
imbalance of transverse momentum from emission of an indefinite number of soft gluons. 
We proposed 

ABNib) =N j rf2p(K,)e*^-^ (13) 

with the condition J (PhAsNib) = 1- 

The soft gluon transverse momentum distribution can be written as 

d^P{\^t) = rf'Ki j .^exp[-l\^t ■ b - h{h, Qma.)] (14) 

where qmax is the appropriate scale for single soft gluon emission and the regularized function 
h{h, qmax) describes the overall soft gluon spectrum integrated over all possible values, up to 

Qmaxy 1-6. 

Hb, qmax) = J - e"^*-^] oc d'k,^[l - e^^*-^] In ^ (15) 
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Eq. (14) was proposed in [37] for soft emission. In its application to hard QCD processes, 



such as lepton pair production in [3H1 EH], an expression for h{b,qmax) was proposed in 



which the integral of Eq. (15) extends only down to a QCD scale /i. This allowed the use 
of the usual Asymptotic Freedom expression (AF) for ag (oif^), safely neglecting the second 



term in the square bracket in Eq. (15). The proposal in [3H] has been used extensively in 
phenomenological analyses and can be traced back also to |10], and the resulting expression 
for (PP(Kt) is usually referred to as the Sudakov form factor. For processes dominated by 
large distances, we have proposed that the kt integration be extended down to zero |11J. Of 
course this calls for the use of an appropriate form for the coupling of soft gluons to the 



emitting quarks, when fcj — )■ in Eq. (15). This strong coupling in the infrared region must 
be such as to give rise to an integrable spectrum. We have discussed this in our papers, in 
particular in [101 [I2] and more recently in |42j . 



From Eqs. (13) and (14), our expression for A{b, s), the normalized Fourier transform of 



the soft gluon resummed transverse momentum distribution, is 

^i?7v(&,^)= ^^,^^.,(,^^_) (16) 

where the dependence of AsNib, s) on the center of mass energy arises through the kinematic 
quantity qmax- In our simplified version of the model proposed in [10], this is obtained by 
averaging over the parton densities. Thus, the rising -high energy part- of the average 
number of collisions, is written as 

nhard{b, S) = ABN^b, S)ajet{s,ptmin) (17) 



B. Scales and parameters 

While we refrain from giving many details of the model, as they have been discussed 
carefully in our different publications [HI [121 [2H] , it is necessary that we give a brief discus- 
sion of the crucial parameters and the different scales involved, so that the physics issues 
become clear. Along with pt, the transverse momentum of the partons involved in the basic 
2 — > 2 scattering, another transverse momentum variable of relevance is kt, the transverse 
momentum of the soft gluons emitted during the hard scattering process. This somewhat 
artificial separation, allows a clear delineation of different theoretical issues involved. The 
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role of ptmin in the calculation of the hard part of the partonic cross-sections has been al- 
ready explained. For the variable k^, an important role is played by Qmax, the maximum 
value allowed to single soft gluon emission during the collision between two partons. Using 
the kinematics of single gluon emission, this quantity can be written as 



qmaAs,y,Q') = ^(1 - — ) , ^ (18) 



2 s a/1 + z sinh^ y 

Here, is the squared invariant mass of the outgoing parton pair, each parton with trans- 
verse momentum pt. In the no-recoil approximation, z = Q'^/s and y is the rapidity of the 
outgoing partons [43j. This scale, which affects the final result mostly through logarithms, 
is a semi-hard scale akin to ptmin and it is also of the same order of magnitude. 

Depending on the energy scale of the process, as discussed at length in |TT] , the calculation 
of the emission of soft gluons of momentum kf < Qmax requires a special treatment. In our 
model, in order to be able to describe very large distance contributions to the total cross- 



sections, we have retained the second term in the square bracket in Eq. (15). Thus, yet 
another scale, in the problem at hand, is the impact parameter b and special care is needed 
in discussing the effects in the infrared region when ktb ~ 1. When b becomes very large and 
kt < l/b < Aqcd, the calculation of h{b,qmax), and hence the overlap function, requires an 
ansatz for the behaviour of in the far infrared where one can not use the usual expression 
at large scales. In the infrared region, our ansatz is to use 

«.// = «f = <^(^)'' (19) 

namely an expression which is singular but integrable, provided p < 1. We stress here that 
to access very large distances, we consider the very soft gluons for which kt — 0, which in 
turns implies resummation and, in the continuum limit, integration. Thus what matters is 
not the limit of the coupling for one soft gluon, but the integrability, a concept also put out 
by Dokshitzer [H] and used in jet analysis 



As discussed below, the expression in Eq. (19) mimics the confinement dynamics. To 



perform the full integral of Eq. (15), which spans from the infrared to the beginning of the 



AF region, we then use an interpolating expression for as, which reduces to the correct AF 



limit at large scales and to Eq. (19) in the infrared region |10j . 

The parameter p captures the physics at scales kt < ^qcd and affects the high energy 
behaviour of the total cross-section in a complex manner. As shown in [12] , its value has to 
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be between 1/2 and 1. As a consequence, the very large b-limit of the function ABNib,s), 
i.e. 

~e-'^('''^) -e-^^')'" (20) 

shows the impact parameter distribution falhng at its fastest as a gaussian {p = 1) and 
at its slowest as an exponential (p = 1/2). The scale A oc Aqcd includes a mild energy 
dependence through the scale qmax, as well as a residual dependence upon the parameter 
p. This behavior in impact parameter space joined with the high energy behaviour of the 
mini-jet cross-sections, ~ s^, was shown to lead to an asymptotic behaviour of the total 
cross-section consistent with the Froissart bound, namely 

^T-|^[eln^]V^ (21) 

upto leading terms in (Ins) [12] . 



1. About the singularity parameter p 

Let us briefly discuss the physics implications of the choice of a particular value of the 
singularity parameter p. In our model, p describes how singular cte// is in the IR region and 
can be related to confinement dynamics. This is seen through the spatial-potential obtained 
through the Fourier transform of the one-gluon exchange potential generated by our effective 
coupling, namely 

V{r) ~ r^P-^ r^oo (22) 

For p = 1/2 the dressed one-gluon exchange potential is (essentially) a constant , whereas 
for p = 1 it is linearly rising. In [l2], we showed that the parameter p has to be between 1/2 
and 1, so that the corresponding potential is confining and Oe// at small scales is singular 
but integrable. On the other hand, given our ignorance of the actual confinement dynamics, 
we shall use it as a free parameter which can interpolate between a fully confining potential, 
linearly rising when p = 1, and the inelastic scenario, in which partons are free to hadronize, 
with p = 0.5 and the potential constant at infinity. 

What we have given above is certainly too simple an argument, but it is not in con- 
tradiction with our understanding of the role played by the infrared singularity. From our 
phenomenology of the total cross-section [261429] . we have found that the value p ~ 3/4 gives 
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acceptable descriptions of high energy data. Such a value would then be consistent with a 
rising potential. 

2. About the PDFs 

We consider the PDFs in our model to be phenomenological QCD parameters through 
which we describe the two high energy effects which are input to the eikonal: the rise at 
the beginning and the slowing down towards asymptotia. So far, ours is a Leading Order 
(LO) description of these two effects. In our parametrization of these effects we use LO 
densities able to go down to very small x- values and ~ 1 GeV^ (see below). In our first 
application of the model in [llj, we have found that GRV densities |16], when combined with 
our resummation expressions, provide a parametrization of the mini-jet effects which leads to 
a good description of total cross-section data. Through the years, other LO PDFs and many 
NLO and NNLO pametrizations have appeared. To use some of them in our model would 
be inconsistent and imply major modifications, as we explain below. Yet another possibility 
will be to employ the LO* distributions |17], suggested for use with the LO Monte Carlos. 
In future work, this and other possibilities will be examined. Here we present the reasonably 
good LO description of the total cross-section, obtained so far in our model in terms of the 
various inputs such as LO PDFs, Ptmin and p, and show how far it can take us in probing 
the long distance behaviour of QCD. 

To be of significance in QCD descriptions of scattering processes, mini-jet models should 
use Parton Density Functions determined from Deep Inelastic Scattering. There are many 
PDFs, and not all of the available sets can be used in a model like ours, in which the rise 
of the total cross-section is due only to mini-jets. To have the rise appearing already at 
ISR energies, one needs Ql = p^min — ^ GeV^ and the PDFs we use must be valid down to 
such a low Qq value. Notice that we can only use LO PDFs, because in our model, the bulk 
of NLO contribution comes from soft gluon resummation, and using NLO versions might 
produce a double counting. The chosen PDFs must also be able to describe very low values 
of X, the energy fraction carried by the incoming partons. In particular most important 
are the low-x gluon densities. Notice however that as the energy increases and lower and 
lower X- values are accessed, soft gluon resummation softens and tames the rise of low-x gluon 
contribution. This contribution is not included in the low-x phenomenology of current PDFs 
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and this may pose a problem since the low x-behaviour of existing densities is modified by 
our model resummation inputs. The claim for instance that the low x-behaviour of some 
of the densities such as GRV is wrong, needs to be seen in light of resummed contributions 
like the ones we are considering : soft gluon effects are a NLO effect, but resummation to 
all orders down into the infrared is not included in the low-x behaviour of PDFs, neither in 
the LO nor in the NNLO parametrizations. Since such effects, in our opinion, are crucial for 
inclusion of QCD in descriptions of the total cross-section, we have resorted to use only LO 
parametrizations for all QCD effects. Even resummation is treated to LO order, since the 



argument of the integral in Eq. (15) does not include NLO terms and Qmax is also averaged 
only over leading valence quak effects (see below). On the other hand, as discussed in the 
Introduction, the bulk of non-leading effects in the resummation is incorporated in our non- 
perturbative ansatz for the singular, but integrable, coupling of soft gluons in the infrared 



region, Eq. (19). 



Our model uses the same approach also in dealing with the calculation of the quantity 



Qmax- Using the values y = 0, pt = Ptmin, "we average the quantity given in Eq. (18) over the 
parton densities, i.e. we use 

,„ , . _ v/i s.,jft^/./.(^i)jtf/i/^fa)v5^j:,..„rf--(i-.) 

2 E.jf^/./«(:^.)/t/,/^(x,)j:L(<i^) 

with Zmin = ^vltminl^- discussed in |12], we only include the leading soft gluon resum- 
mation from valence quarks and not from low-x gluons. The emission of soft radiation only 
from valence quarks, and thus the determination of Qmaxi is also part of our LO parametriza- 
tion, and follows the idea that in resummation in QED, the leading terms correspond to the 
external legs. Although in the impulse approximation gluons are treated as free particles, 
gluons in gluon-gluon collisions must partake of the initial acollinearity imparted by the 
valence quarks. It would be inconsistent to imagine that the acollinearity due to initial state 
radiation from valence quarks does not reflect on all subprocesses. This is why at LO we 
consider this to be the major effect. The introduction of emission from interacting gluons 
will have to be considered but the formalism is more complicated, and we flrst want to see 
how far the LO model can take us. 

Among the available PDFs satisfying these requirements, there are GRV, MRST and 
CTEQ densities. The latter, in our model, do not yield a rising cross-section past the 
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TeVatron data. Within our approach, the reason hes in the fact that the quantity qmax, 
calculated with CTEQ, rises with energy so much that the taming effect from soft gluons, 
which reflects the initial acollinearity of the partons, does overcome the rise of the mini-jet 
cross-sections. We have shown the details of this exercise in |18] and work is in progress 
towards a better understanding of how to use CTEQ in our model. To summarize, our 
choice of PDFs is a parametrization of LO QCD effects within an eikonal mini-jet model 
with soft gluon resummation in the infrared region: as such, densities like CTEQ do not 
describe the total cross-sections and cannot be used in the model in its present formulation. 
On the other hand, both GRV as well as MRST LO densities are well suited to describe all 
the data up to the Tevatron and including cosmic rays, and can thus be used for further 
investigation. 

NLO corrections to our model would most certainly change the choice of model param- 
eters, of course. An NLO version of our model might clarify some of the issues and a full 
investigation in this direction is in progress. 

Thus, at the end, the relevant parameters for evaluation of fihardip^s) are ptmin,P and 
the chosen PDF set. We call them the high energy (HE) parameters. As noted, the scale 
Qmax, the upper limit in the soft kt integration, is not an independent parameter. Its value is 
determined by the kinematic variables defining the parton-parton sub-process, namely pt, the 
rapidity y and s. In our model, for the calculation of Qmax, we choose y = 0, pt = Ptmin- The 
scale Qmax, uscd in obtaining the impact parameter distribution ABN{b,s), is then calculated 
as an average over the parton densities in the colliding hadrons. Once the parameters ptmin 
and the PDF set are chosen, ajet and Qmax can be calculated. One then chooses a value for 
the parameter p, and using qmax, one can completely determine ABN{b,s). The calculation 
of nhard{b, s) then follows from choosing p, ptmin and the PDFs [101 flT] - 

IV. THE TOTAL CROSS-SECTION 

The mini-jet model we have outlined in the previous section describes only the rising part 
of the total cross-section. While the hard term ntardib, s) dominates at high energy and is 
supposed to include the dynamical mechanism of the rise of the cross-section with energy, 
a substantial part of the total cross-section is also due to soft processes with pt < ptmin, 
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according to our model. Due to our limited understanding of the QCD dynamics of this 
part of the cross-section, nsoft{b,s) is parametrized. We start with factorization 

nsoft{b, s) = Asoft{b, s)asoft{s) (24) 

and have used two models to calculate Usoft- The two models use two different expressions 
for the overlap functions and different ansatze for the energy dependence of the soft cross- 
sections. 

The details of these two models are discussed in Appendix [A], where we also give values 
of the parameter sets we use. We stress here the important fact that our parametrizations 
for Usoftib, s) do not include any term rising with energy, but only constants and decreasing 
terms. Nowhere we impose in the model any logarithmic growth. The rise is always produced 
by the mini-jets. 



Notice that Eq. (10) can be written as 

n(6, S) = Amean{b, s)[cr^o/t(s) + (^jet{s)] (25) 

with the overall mean impact parameter distribution defined as 

^mean{0,S — ;r^— -- {AO) 

CTsofAs) + (Tjet[s) 

thus allowing for comparison with other models. In the four plots of Fig. [T] we show the 
model prediction for this function and other impact parameter distributions of relevance 
to the total and inelastic cross-section calculation. In all four figures, we used Model II for 
^softib, s). Using Model I, gives very similar results. 

In the two upper plots of Fig. [l] we show the behaviour of Ameanib, s) and of the aver- 
age number of collisions, n{b, s) for proton processes as a function of b at ^/s = 7 TeV. 
Comparison is shown between the overall distributions and the contribution from soft and 
hard components of the imaginary part of the eikonal function. The bottom plots show the 
elastic amplitude in the approximation ^ex{b, s) ~ 0, and the integrand for the inelastic 
cross-section calculation in the two dimensional b space. These plots (bottom) clearly show 
the cut-off behaviour discussed in [28], resulting from our soft /Ef-resummation model in the 
infrared region. It is interesting to compare, in the bottom left figure, the behaviour at low 
cm. energies with the high energy LHC behaviour, which is more and more approaching a 
Fermi function in 6-space.This cut-off behaviour is what transforms the power-like behaviour 
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b (GeV"'') b (GeV"'') 

FIG. 1: The impact parameter dependence of different quantities relevant for total cross-section 
calculations are shown. Top left : the normalized Fourier transform of the soft gluon /ct-distribution, 
which gives the overlap function in proton-proton collisions, averaged over densities, as well as over 
hard and soft processes. Top right: the average number of collisions as a function of the impact 
parameter. Bottom: at left, elastic scattering amplitude in b space, in the approximation Ke^ ~ 
at various cm. energies, and, at right, the probability for inelastic collisions at y/s = 7 TeV, 
namely the integrand for the inelastic cross-section. The high energy parameters are ptmin = 1.15 
GeV, GRV densities I46j and p = 0.75, except for the dashed curve in the bottom right panel, 
where p = 0.5. 

of the mini-jet cross-sections into a logarithmic energy dependence coherent with the limits 
imposed by the Froissart bound |12] . 

The bottom figure at right illustrates the role played by the parameter p in our model. 
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We notice that the plateau extends to include larger 6— values as p decreases from the value 
p = 0.75 to p = 0.5, indicating a larger interaction region when p is smaller, i.e. when 
the infrared singularity, accessed by soft gluon resummation, lessens. In other words, as p 
decreases from 0.75 to 0.5, the less confining behaviour for soft gluons in the infrared enlarges 
the interaction region, whose extension is roughly determined by the value where the plateau 
drops. In Sect. |V], we shall compare the data for amei = o-totai — o'eiasUc with the inelastic 
cross-section from the model, Eq. Q, and we shall find that the value p ~ 0.75, which 
describes total cross-section data, falls short of the totality of inelastic data for energies up 
to the TeVatron. We shall then resort to a smaller p value to describe them. 

We close this discussion of the distributions in impact parameter space, with a few 
general remarks. As is evident from the bottom curves in Fig.jlJ at high energy the integrand 
in impact parameter space is a constant until a certain value bmax{s) (approximately about 
one Fermi at ^/s = 7 TeV), and then there is a tail from the peripheral collisions. As the 
initial energy of the two protons increases, the plateau is extended and the tail becomes 
sharper. In the extreme high energy limit, the distribution is expected to become sharp (a 
perfect Fermi function) and in this perfect hard disk limit the elastic cross-section should 
become one-half of the total cross-section. While the trend of the data so far is towards an 
increase in this ratio, the approach is extremely slow and not likely to be reached in the 
near future. 

We are now in a position to exhibit our results for the total cross-section. In Fig. |2] we 
show two black lines and a band: the full black and the green band are obtained with Model 
I for Usoft and have been already presented [III [27] . The dashed black line was obtained 
with Model II, and the low energy fit includes the mini-jet contribution. For all the fits, 
we have used LO PDFs which allow values as low as can be reached for the chosen 
Ptmin well as the corresponding low x values that are reached. As discussed earlier, we 
use here the following parameterisations: GRVUg, GRV94 [49j, GRV98 [SO], MRST 
The two curves, Model-I and Model-II, use the same set of high energy parameters, GRV 
densities for the mini-jets, ptmin = 1-15 GeV and p = 0.75 for the infrared singularity 
parameters. The green band is obtained varying slightly the high energy parameters, by no 
more than a few percent, and using different partonic densities. At present LHC energies, 
i/i = 7 TeV, we obtain atotai = 91.6 + 6.5 — 9.7(m6) with Model I . The results for 
this model at the two LHC energies y/s = 7, 14 TeV, are summarized in Table m For 
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E 




FIG. 2: Total pp and pp cross-section calculated with our model and compared with data. Values 
for atotai at 7 TeV at LHC are indicated. 

Model II, we have only calculated the value with GRV densities, ptmin = 1-15 GeV, and 
p = 3/4 , and the value (Jtotai{y/s = 7 TeV) = 93.5 mb, lies within the indicated band of 
Model I. Our actual ignorance of a fundamental QCD description of the low energy total 
cross-section dynamics makes further studies of the difference between the models, quite 
irrelevant. However the result with low energy fits inclusive of mini-jet contribution, i.e. 
model II, should be preferred. 

Having thus chosen the set of HE parameters which adequately describe the existing total 
cross-section data, we now proceed to investigate how the eikonal function, determined by 
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PDF 


Ptmin (GeV) 


CTo(mfe) 


p 






atot (mb) 












= 7 TeV 


^ = 14 TeV 


GRV 


1.15 


48 


0.75 


91.6 


100.3 


GRV94 


1.10 


46 


0.72 




95.9 


103.9 




1.10 


51 


0.78 




83.3 


89.7 


GRV98 


1.10 


45 


0.70 




94.9 


102.0 




1.10 


50 


0.77 




82.1 


87.8 


MRST(72) 


1.25 


47.5 


0.74 




85.9 


95.9 




1.25 


44 


0.66 




97.9 


110.5 



TABLE I: Values of atot for ptmin, co and p corresponding to different parton densities in the 
proton, for wliich our model [27] gives a satisfactory description of the total cross-section. The 
values are obtained with Model I (see Appendix) for Ugoft^ the average number of soft collisions. 

these parameters, describes data for the inelastic cross-section. 

V. THE INELASTIC PROTON-PROTON CROSS-SECTION 

In this section we present the results of our model for the inelastic cross-section. It 
is worth noting that a nontrivial bound on the inelastic cross-section, based on general 
analyticity arguments has been obtained only very recently ||52j, unlike its counterpart for 
the total cross-section which has been known for quite some time 

The inelastic proton-proton cross-section is of particular interest in cosmic ray physics, 
as it determines the total p — air cross-section [33]. The knowledge of this inelastic cross- 
section at very high energies enters the simulation of cosmic ray air showers and is used for 
extracting the total pp cross-section from cosmic ray data. However, unlike the total and 
the elastic case, the inelastic cross-section is not uniquely defined and data for it suffer from 
usage of different cuts imposed in the analysis which lead to inlcusion of differing amounts 
of diffractive contribution to the measurement. In principle, the least ambiguous way to 
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define data for tlie inelastic total cross-section is through the definition 



exp exp exp (0'7\ 

''^inelastic ''^ total '^elastic \ ) 

Models can then be compared with data in the available range of energies, starting from 
^/s ~ 1 — 2 GeV up to the TeVatron data, and now to the LHC data at ^/s = 7 TeV. 

As a first step towards analysing the inelastic cross-sections, one then needs to extract 
the values for pp scattering, as the difference between the total and the elastic cross-sections. 
The error of this procedure is obtained by combining the errors in quadrature. Once the 
data points are obtained in this fashion, one can try to confront our model predictions with 
them. Note that there are no data for pp scattering beyond the ISR energies [56j . However, 
the increasing dominance of gluon-gluon scattering with the rising energy should allow us to 
use the pp data from CERN SppS [57] and from the TeVatron [55H^ . as a guidance in our 
analysis. We notice that in [58], in addition to the total and elastic data, the collaboration 
also presents a value for the inelastic cross-section at ^/s = 1800 GeV. There are therefore 
two possible values to fit, atotai — (^elastic = 56.2 ± 3.5 mb (taking atot from |59] and o"ez from 
[58] ) or a inelastic = 55.5 ± 2.2 mb which is the direct measure presented in [58]. The two 
values agree within the errors. In our estimate, for consistency with all the other data points, 
which were obtained in the same way, we use the value obtained by subtraction. For the 
same reason, i.e. consistency with the subtraction procedure we have outlined, we do not 
use results from [63j. Thus, in all the figures to follow, the term inelastic data corresponds 



to data obtained from the difference between total and elastic cross-section, as per Eq. (27) 
and can be directly compared with Eq. (|4]). 

As was the case for the total cross-section, before studying the high energy behaviour, 
the soft part of the eikonal needs to be chosen. We shall use the Form Factor model. Model 
II, with two different choices for the soft cross-section asoft- In the first. Model II-A, the 
soft part of the eikonal, namely n^o/t = AppO'soft, is the same as the one entering the total 
cross-section estimates shown by the dashed line in Fig[2| while in the second. Model II-O, 
we use an independent fit to the low energy inelastic pp cross-section data with ^/s > 5 
GeV, which also includes the mini-jet contribution. As in the case of the total cross-section, 
at high energy, the cross-section behaviour is rather insensitive to this choice and is mainly 
controlled by the PDF's and the other high energy parameters p,ptmin- 

As anticipated in Sec. [11} we find that the high energy data, cannot be described by the 
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FIG. 3: The total and inelastic cross-section as described by our model when using the same set of 
high energy parameters in the eikonal. Results are compared with pre-LHC data and an estimate 
from 15 1 . 



same set of HE parameters determined from Eq. ([7|) for the total cross-section, namely 
GRV and MRST PDFs, ptmin and p. In particular, in Fig. [3] we show a comparison 
between the inelastic data and the results of our model for the set of HE parameters GRV, 
Vtmin = 1-15 GeV.p = 0.75. More generally, varying ptmin and the densities as in Table 
[T] within the same ranges as we did for the total cross-section, we find that the p-values 
in Table [T] always give results short of the pre-LHC data. This is consistent with the 
observation in Sec. [n| namely that the eikonal mini-jet model with two components will only 



describe uncorrelated processes, while the inelastic data from Eq. (27) include also collisions 



in diffractive and other similarly correlated regions. On the other hand, everything else 
being equal, a reasonable description is obtained by letting the parameter p vary between 
the minimum value consistent with the model |12], p = 0.5, and 0.66, the highest value 
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Vs ( GeV ) -Js ( GeV ) 

FIG. 4: Predictions for the inelastic pp cross-section from model II-A (left) and Model II-O (right). 

which can possibly be consistent with the total cross-section band in [27]. We show this 
analysis in Fig. |4j 

We combine the results of Fig. |4] in Fig. |5| In this figure, using both Model II-A and 
Model II-O to constrain the low energy behaviour, we indicate all the parameter sets used 
to produce the bands, as well as the numerical values of interest. We see that the low energy 
behaviour is of course better described by Model II-O, which fits low-energy inelastic data, 
including mini-jets in the overall fit. However, a good description is also obtained from the 
same low energy eikonal entering the total cross-section (Model II-A), an indication that 
the two-channel eikonal model works at low energy. In this figure, for simplicity, we show 
the application of Model II-O only for GRV densities. Also, to compare our choice of high 
energy parameters with those used for the total cross-section, we have plotted a curve (black 
line) with same ptmin = 1-15 GeV and GRV densities as in Fig. [sj but with p = 0.5, a value 
which we find to give a good description of the pre-LHC inelastic data. 

VI. ELASTIC, INELASTIC AND DIFFRACTION PROCESSES IN EIKONAL 
MODELS 

In the previous section, we have seen that, up to TeVatron energies, the two channel 
model cannot accomodate both inelastic (diffraction included) and total cross-section data 
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FIG. 5: Comparison between data for atotai — (^elastic with results of our model. Values for the 
parameters used for the band and the curve are indicated. 

without a change of parameters. On the other hand by changing the singularity parameter 
p, in a way consistent with our understanding of the role played by the parameter, the model 
results can span both the total and the inelastic data at energies up to the TeVatron (blue 
band). We shall now combine the results of our eikonal mini-jet model for the total and the 



23 



inelastic cross-section in a single figure and compare the latter with ATLAS recent results [2] 
and with a recent calculation by Block and Halzen |3] . In Fig. [6] the upper (green) band is 
from Fig. |2] and the lower (blue) band from Fig. |5] . Once a set of parameters has been fixed 
from the total cross-section data, the dashed curves indicate what the two-channel eikonal 
mini-jet model would predict, without changing any parameter in going from the total to the 
inelastic. As already mentioned, the dashed curve for the inelastic cross-section does not fit 
the TeVatron data, nor the SppS data, but its comparison with the recently released ATLAS 
measurement or preliminary CMS data in the central region can be very instructive. The 
dashed line (see Fig. [s]), corresponding to the model with the same p value as in the total 
cross-section curves, is close to the preliminary value reported by the CMS Collaboration for 
inelastic events with at least 2 charged particles with t] < 2.4 and p± > 200MeV/c^ [3J This 
confirms the interpretation of the two-channel eikonal minijet model, put forward in Sect. 
II, namely that this model describes total inelastic processes, with little or no correlations. 

On the other hand, with the present parametrization, but changing the parameter p, 
we find that our lower (blue) band is consistent with the ATLAS results, as corrected for 
diffraction |2] (full circle). 

Let us now discuss the missing element in Eq. Q, that is a category of processes which 
are inelastic but correlated, namely single and double diffraction. In diffraction, partons in 
the final state resulting from multiple collisions or from soft gluon emission must recombine 
into a proton or are constrained to be emitted along the outgoing proton. It would then be 
more correct to write Eq. ( [6| as 



^uncorrelated^^^ ^ f d% P{{n,n}) = I £h[l - e'^^'''^)] (28) 

n=l 

SO that, generally speaking, 

_ I diff I uncorrelated fOQ\ 

(y total - (y elastic + O^inel + ^inel ['^^) 

and include in crf^{/ all correlated processes, like single and double diffraction. Thus Eq. Q 
can be used only for inelastic uncorrelated processes, while the inelastic correlated processes 
are to be described separately. 

On the other hand, the above considerations do not yet suggest how to include the 
diffractive processes. The solution adopted in this paper, illustrated by the lower (blue) 
band in Fig. [6} has been to mimick the presence of diffraction processes using different 
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FIG. 6: Total and inelastic data for pp and pp scattering compared with the prediction from our 
model (dashed curve) [3], in the two-channel case, when the same eikonal function is used both for 
the total and the inelastic cross-sections, from low to high energy. The upper band is the same as 
in Figj2j the lower band for the inelastic cross-section is the same as in the left panel of Fig. [4} 
Comparison is made with theoretical predictions by Block and Halzen [5], with TOTEM data [8] 
for the total cross-section, ATLAS inelastic data extrapolated to ^ > rrip/s and with ATLAS data 
in the range 5 x 10^^, where ^ = M'j^/s (star symbol) [2]. 

values for the singularity parameter p in Eqs. ^ and ([T]). By comparing with previous 
data, Figs. [5] and [6] indicate that the value p = 0.5, with the present LO parameterization 
(LO GRV and MRST), gives the best description of the total inelastic cross-section, while 
the above discussion, following the one in in Sec. |TT], indicates that values with p = 0.75 
would only include uncorrelated inelastic processes. 

For the future, once the parameters of the model have been tuned to the measured total 
cross-section, more precise predictions for the totality of uncorrelated inelastic events could 
be extracted from our model. 
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A. Uncertainties and the nature of the theoretical band 

For clarity and completeness, we recapitulate here the uncertainties in our model and 
the nature of the bands we have presented. The bands in this paper preceded the release 
of LHC data for total and inelastic cross-sections, and reflect how well the model worked in 
describing lower energy data [H |3T] as well as their uncertainties, mostly the TeVatron data 
for the total cross-section. For the future, the new data will allow a better tuning of all the 
parameters involved and hence firmer predictions for total cross-section data at ultra high 
energies or for quantities, such as survival probabilities, needing the inelastc cross-section. 

Even so, the bands we have presented are in reasonable agreement with the new data 
and provide confirmation of our understanding of the eikonal two-channel model. The blue 
(lower ) band in Fig. |6] describes the interplay between all the HE energy parameters. The 
essential parameters are the PDFs, ptmin and p. The PDFs control both the rise (mini-jets) 
as well as the logarithmic asymptotic behavior through Qmax- At fixed Ptmini different PDF's 
lead to different values for a jet and for qmax- In the present formulation, a jet are controlled 
by uncertainties in the LO low-x behaviour of the gluon densities, and q^ax by the LO 
valence quark densities. The lower is ptmin the faster the rise, and, at high energies, the 
more this rise needs to be quenched. The quenching of the rise depends on how much is 
the allowed loss in transverse momentum, qmax and upon the degree of singularity of the 
effective coupling, the parameter p, for which we have both an upper as well as a lower 
bound. For the convergence of the IR integrals, it is incumbent that p < 1. [It must be 
recalled that the Wilson equal area rule leading to a strictly linear potential viz., p = 1 is 
not allowed in our resummed model] . On the other hand, it is easy to check that there is no 
confining potential if p < 1/2. 

Within the contraint 1/2 < p < 1, values for this parameter are determined through 
the phenomenological description of the data. There is an interesting interplay between the 
phenomenological value of p and its weight in describing the uncorrelated versus correlated 
events. In the model, once the value for p has been fixed through the total cross-section, it 
can be used to describes the uncorrelated inelastic events, while the fully inelastic data need 
a different treatment. We have mimicked the amount of diffraction to include in the full 
inelastic cross-section by decreasing values of p. As p is decreased, the amount of acollinearity 
introduced by the soft-gluon radiation changes. We see, that a description of the underlying 
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fully inelastic data would favour p = 0.5. 

Finally, we point out that some uncertainties are also introduced by the different purely 
phenomenological models for the quantity fisoft and are reflected in the band for the inelastic 
cross-section. 

VII. CONCLUSIONS 

We have presented results from the two channel eikonal model for the inelastic cross- 
section. We have shown that such formulation only describes uncorrelated processes, a 
result independent of the details of the model used to describe the data. This observation 
clarifies the origin of a long standing difficulty of the two channel eikonal models well known 
to people working in total and inelastic cross-sections. 

We have shown numerical results for the description of total and inelastic cross-sections 
at LHC at ^/s = 7 TeV and beyond, using our QCD based model. Our numerical results are 
obtained in a LO QCD model for collisions, embedded in the eikonal formulation. The major 
QCD phenomenological input is from Parton Density Functions, which can be considered 
QCD parameters through which we build our eikonal in order to describe the two high 
energy effects: the rise at the beginning and the slowing down towards asymptotia. While 
it is difficult to describe these two effects in QCD models which satisfy unitarity and the 
Froissart bound, in our model we have both. 

We have found that with exactly the same low energy parameters entering the total cross- 
section calculation, but with different values of the high energy parameters p and ptmin, our 
model predictions catch the TeVatron data as well as new data from the ATLAS collabora- 
tion [2j . On the other hand, if no parameters are changed at all in the eikonal function used 
for (Jtotah our results agree with the hypothesis that the simple two-channel eikonal mini-jet 
models describe only the non-diffractive part of the inelastic cross-section. We defer to fu- 
ture work an analysis of the inelastic cross-section which might include correlations in the 
diffractive terms, or inclusion of a non zero 3fJex(&, s). 



27 



Acknowledgements 

We acknowledge conversations and discussions with L. Anchordoqui and thank him for 
providing us with the numbers we have used for the inelastic cross-sections. G.P. thanks 
the MIT Center for Theoretical Physics and Brown University Physics Department for hos- 
pitality while this work was being written. Y.N.S. thanks Northeastern University Physics 
Department for hospitality. Work partially supported by Spanish MEC (FPA2006-05294, 
FPA2010-16696 and ACI2009-1055) and by Junta de Andalucia (FQM 101). This work has 
been supported in part by MEC (Spain) under Grant FPA2007-60323 and by the Span- 
ish Consolider Ingenio 2010 Programme CPAN (CSD2007-00042). R.M.G. wishes to ac- 
knowledge support from the Department of Science and Technology, India under Grant No. 
SR/S2/JCB-64/2007, under the J.C. Bose Fellowship scheme. She also wishes to thank the 
theory division at CERN during the preparation of an earlier version of this manuscript. 



Appendix A: The average number of soft collisions nsoft{b,s) 

As in the case of the hard collisions, we approximate nsoft{b,s) with the factorized ex- 
pression 

nsoftib, s) = Asoft{b, s)asoft{s) (Al) 

and then use two different phenomenological parametrizations for the average number of 
soft collisions. In the first model. Model I, we assume that the decrease in the proton proton 
total cross-section is due only to soft gluon emission. 



a. Model I, Asoftib,s) given by resummation expression: 

This is the model adopted in our analyses in Refs. [121 EZ]- In this case, for the b- 



distribution, we use Eq. (16), i.e. the BN model. In this equation, we parametrize qmax 
to make it vary slowly with energy, thus giving a slight (decreasing) energy dependence to 
Asoft{b,s) [12]. In Table [n] from [12] we reproduce the values of g^/* used to evaluate 
Asoft{b,s). These values correspond to the observation that for processes contributing to 
Usoft, a soft gluon will always carry away less energy than for those contributing to Uhard- 
The question is how much lower is the allowed energy. The above set of values is purely 
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TABLE II: Average qmax values used for the impact parameter distribution of the soft part of the 
eikonal 





qmli{GeV) 


5. 


0.19 


6. 


0.21 


7. 


0.22 


8. 


0.23 


9. 


0.235 


10. 


0.24 


50. 


0.24 


100 


0.24 



phenomenological and the values are those we found to give an acceptable description for 
(Tpp before the rise. 

For pp scattering, where the s-channel is exotic and no trajectories are exchanged in the 



t— channel, our use of soft gluon emission using a form for AsNi^^ s) a-la Eq. (16) with q^ax 
chosen as explained above, is sufficient to explain the slight decrease at low energies, in the 
total pp cross-section. The situation is more complicated for pp. As described in p^, we 
then simply write for the cross-section 

2e , 



with (To a constant, as in Table |T]and e = 0, 1 for pp and pp respectively. 



(A2) 



h. Model II, Asoftib,s) given by ApF{b,s): 

In the second approach, for the soft collisions we have used the Form Factor model for 



A(b,s), Eq.( 12), together with a standard, Regge exchange inspired parametrization oi pp 



and pp cross-sections at low energy, namely 

^""W=''°+(£i^=Fp;^ (A3) 

with =p referring to pp and pp respectively, and the constant parameters determined through 
a fit to the low energy data. In earlier publications, such as Ref. [11], the fit to low energy 
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TABLE III: Results of the fit to pp and pp data. 



Fit for pp 


Fit for pp 


ao = (48.20 ±0.19) mb 
^1 = 101.66 ±16.35 ai = 0.99 ±0.13 
^2 = 27.89 ± 4.78 02 = 0.59 ±0.04 
= 154.1/(102±5- 1) 


ao = (47.86 ± 2.47) mb 
yli = 132.07 ±32.89 ai = 0.69 ±0.14 
yl2 = 0.82 ± 0.31 02 = 0.52 ±0.07 
= 24.65/(31 + 5 - 1) 



data was done using an expression which did not include the mini-jet contribution, whereas 
in our more recent analysis |28] the minijet contribution is included in the fit. From the 



point of view of the x , the difference in the results is not noticeable. In Table III from [28] 



we reproduce the values of various parameters used in Eq. (A3) 
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